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ABSTRACT 


This  paper  evaluates  five  maneuvering  strategies  for  am 
evader  in  a  two  dimensional  continuous  evasive  game.  An 
evader’s  movement,  in  two  dimensions  at  constant  speed,  is 
simulated  by  choosing  courses  that  are  independent  and  dis¬ 
tributed  according  to  one  of  the  five  course  change  rules „ 
The  times  between  course  chances  are  independent  exponential 
random  variables.  An  improvement  over  a  previously  estab¬ 
lished  least  upper  bound  on  evader  survivability  is  pre¬ 
sented.  The  bulk  of  this  paper  is  the  presentation  of 
evader  survivability  as  a  result  of  the  maneuvering  strategy 
he  employs. 
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I.  INTRODUCTION 


A  classic  military  problem  with  current  application  is t 
how  to  maneuver  a  mobile  target  in  order  to  prevent  success¬ 
ful  prediction  of  its  position.  For  example  how  should  s  i 
aircraft  carrier  maneuver  to  maximize  its  survivability 
against  an  ICBM  attack? 

Of  primary  importance  to  this  problem  is  the  time  lag 
between  the  decision  to  fire  at  a  target,  made  at  time  t, 
and  the  warhead  detonation  at  time  t+TL.  This  lag  is  called 
time-late  and  denoted  as  TL.  The  decision  to  initiate  such 
an  action  is  based  in  part  upon  knowing  the  target's  posi¬ 
tion  at  time  t  while  the  effectiveness  of  that  action  is 
dependent  upon  the  error  in  predicting  that  position  at  t+TL. 

A.  ICBM  VERSUS  AIRCRAFT  CARRIER 

Consider  the  ICMB  versus  aircraft  carrier  problem.  As¬ 
sume  that  the  carrier  operates  in  raid-ocean,  possesses  great 
endurance,  can  make  sharp  turns  and  has  as  its  only  kinematic 
restriction  a  constant  speed,  v.  The  carrier  is  referred  to 
as  the  evader,  E.  The  attacker,  P,  continuously  observes  E 
from  a  nearby  unarmed  trawler.  P's  weapon  is  a  land  based 
ICBM  which  has  perfect  accuracy,  produces  a  lethal  area  A 
(this  lethal  area  can  have  any  shape)  but  does  not  have  a 
mid-course  guidance  nor  homing  capability.  Although  P  knows 
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E's  position  past  and  present,  he  also  knows  that  his  ICBM 
will  have  a  one  hour  time-late. ^ 

P's  problem  is  to  predict  E's  position  at  time.  t+TL  in 
order  to  initiate  an  attack  at  time  t.  E's  problem,  since 
he  will  not  know  of  an  attack  until  he  observes  the  deto¬ 
nation,  is  to  maneuver  in  such  a  way  to  confound  p's  predic- 
tion.  E's  movement  should  be  random  because  any  nonrandom 
movement  would  be  vulnerable  to  extrapolation. 

B.  AN  UPPER  BOUND  ON  EVADER  SURVIVABILITY 

A  lower  bound  on  P's  success,  p^,  and  consequently  an 

upper  bound  on  E's  survivability,  1-p^,  is  known.  P  knows 

E's  position  at  t+TL  will  be  interior  to  or  on  a  circle  of 

radius  v*TL  centered  on  E's  position  at  t  (which  is  known). 

The  enemy  can  then  gaurantee  a  kill  probability  of  at  least 
2 

A/ir(v«TL)  by  choosing  a  lethal  area  randomly  within  the 
circle  of  uncertainty  (a  wedge  of  random  orientation  will 
do)  [Ref.  4] .  Then  the  maximum  survivability  E  could  attain 
would  be  1-A/tr  (v-TL)  2. 

This  upper  bound  on  E's  survivability  is  a  function  only 
of  P's  strength,  that  is  the  magnitude  of  A  (the  value 
tt(v*TL)  is  constant).  Figure  1  is  a  graph  of  that  survival 
function  versus  P's  probability  of  success. 

^•Time-late  is  a  summand  of  many  factors  in  the  command  and 
control  problem  as  well  as  the  missle  flight  time.  For  the 
purpose  used  here  it  is  sufficient  to  treat  time- late  in 
total  and  not  consider  its  decomposition. 

2 

If  E  knows  the  time  at  which  an  attack  is  initiated  then 
the  problem  is  trivial. 
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Versus  Attac! 


FIGURE  1 
ipecific  set  ■ 


II.  BACKGROUND 


A.  DISCRETE  EVASIVE  GAMES 

It  was  felt  that  a  two  dimensional  continuous  game  was 
too  difficult  to  solve.  So  the  approach  was  to  solve  a 
similar  but  easier  game.  It  was  assumed  that  the  ocean  was 
a  linear  set  of  discrete  points  and  that  E's  mobility  con¬ 
sisted  of  being  able  to  jump  either  right  or  left  to  an 
adjacent  point  [Refs.  2  and  3].  Time-late  became  an  integer 
number  of  jumps  that  E  could  make  between  the  attack  de¬ 
cision  and  warhead  detonation.  As  in  the  exact  problem  E 
had  no  knowledge  of  the  decision  to  attack. 

Discrete  problems  of  this  type  are  called  evasive  games 
and  are  classified  by  the  number  of  jumps  constituting  time- 
late.  A  one  step  discrete  evasive  game  means  E  can  move 
either  left  or  right  one  jump  prior  to  detonation.  Game 
theory  provides  an  immediate  solution  to  that  game;  the 
value  is  p^O.5.  Also  there  is  an  optimal  strategy  for  E 
in  that  game.  At  each  jump  E  should  go  left  with  probability 
one-half  or  right  with  probability  one-half.  Employment  of 
such  a  rule  for  each  jump  guarantees  E  a  survival  probability 
of  at  least  the  game  value,  regardless  of  P's  firing  rule. 

Having  easily  solved  the  one  step  game  interest  was 
focused  on  the  two  step  game.  In  this  game  E  was  allowed 
two  jumps  during  time-late.  The  solution  has  been  obtained, 
but  it  was  not  as  easy  to  achieve  as  that  of  the  one  step 
game  [Ref. 2] . 
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The  direction  of  the  analysis  was  clear.  Knowledge 
gained  from  solving  the  more  simple  games  would  be  a  step¬ 
ping  stone  to  the  solution  of  the  more  difficult  games. 
Eventually  the  assumptions  of  linearity  and  discreteness 
could  be  relaxed.  Progress,  however,  in  solving  the  dis¬ 
crete  games  has  been  very  slow.  Researchers  are  presently 
embroiled  in  solving  the  three  step  game. 

B.  A  CONTINUOUS  EVASIVE  GAME 

Washburn  in  Ref.  1  presented  a  different  approach  to 
analyzing  the  ICBM  versus  aircraft  carrier  problem.  He  de¬ 
veloped  the  probability  density  function  of  a  particle 
moving  continuously  in  a  two  dimensional  medium  subject  to 
a  specific  set  of  maneuvering  rules.  That  approach  was 
distinct  from  previous  work  because  it  addressed  the  exact 
problem.  The  results  presented  in  this  paper  are  an  ex¬ 
tension  of  that  approach. 

Washburn  proposed  the  following  strategy  for  E  and 
analyzed  the  evader's  subsequent  survivability,  first  as  a 
function  of  the  attacker's  strength  and  secondly  as  a 
function  of  a  specific  strategy  parameter.  E  was  to  select 
courses  from  a  uniform  distribution.  Course  changes  were 
to  occur  as  a  Poisson  process.  This  meant  the  time  between 
course  changes  would  be  exponentially  distributed. 

The  exponential  time  between  course  changes  has  an 
intuitive  appeal  because  of  the  memoryless  property  of  that 
distribution.  The  memoryless  property  is:  the  probability 
that  E  will  not  change  course  in  the  future,  given  that  he 
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has  not  changed  course  for  some  observation  period,  is  in¬ 
dependent  of  the  length  of  that  period.  Therefore,  use  of 
the  exponential  time  distribution  should  serve  to  confound  a 
prediction  of  future  position  by  extrapolation,  regardless 
of  the  course  distribution  used.  For  example,  suppose  P's 
attack  decision  rule  required  that  he  observe  E  maintain  a 
constant  course  for  at  least  five  hours  prior  to  initiating 
an  attack.  Such  a  procedure  would  not  improve  P's  prob¬ 
ability  of  a  kill  because  the  only  information  of  any  benefit 
at  the  time  an  attack  is  initiated  is  E's  position  and  last 
course. 

Figure  2  is  a  graph  of  the  evader's  maximum  attainable 

survivability  when  using  the  strategy  of  uniform  courses 

and  exponential  times  versus  P's  strength  for  the  optimal 

3 

exponential  parameter,  X.  For  comparison  the  upper  bound 
on  survivability  is  also  presented. 

"If  we  assume  the  parameters  A,v,TL  are  known  to  both 
sides,  then  the  evader  can  select  X  to  maximize  the  sur¬ 
vivability.  The  evader  will  clearly  be  in  trouble  if  he 
makes  X  too  small,  because  the  kill  probability  is  at  least 
exp(-X»TL).^  On  the  other  hand  making  course  changes  too 
frequently  will  lead  to  a  density  function  that  is  highly 
peaked  at  the  origin,  which  is  equally  undersirable . . . " 

[Ref.  1]. 

^The  exponential  parameter,  X,  is  the  inverse  of  the  mean 
time  between  course  changes . 

^Exp(-X»TL)  is  the  probability  that  E  will  make  no  course 
change  during  TL. 
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Minimum  p^  for  Optimal  Parameter  \  Versus  Attacker 
Strength;  Strategy  Is  Uniform  Courses 
and  exp(X)  Times 

PI GURU  2 


"It  is  not  known  whether  or  not  the  optimal  strategy  for 
E  is  a  Poisson  strategy  of  the  type1  just  considered,  or  even 
whether  the  uniform  distribution  on  angles  is  optimal  within 
the  class."  [Pef.  2] 
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III .  RESULTS  OF  SIMULATING  THE  CONTINUOUS  EVASIVE  GAME 

The  goal  of  this  research  was  to  determine  if  the  uni¬ 
form  course  rule  was  optimal  within  the  class  of  Poisson 

i 

.  ■  strategies  or  to  find  a  better  rule  if  it  was  not.  For 

another  rule  to  be  better  than  the  uniform  one  it  would 

\ 

have  to  enable  E  to  attain  a  higher  survivability  than  that 
shown,  in  Fig.  2. 

The  methodology  of  solving  discrete  games  had  not  ad¬ 
vanced  sufficiently  to  achieve  those  goals.  The  methods 
used  by  Washburn  were  also  not  suitable  because  of  the 
mathematical  difficulty  of  the  problem.  'i 

|  The  Poisson  class  of  strategies  is,  however,  uniquely 

\ 

suited  for  analysis  by  computer  simulation.  This  is  due  to 
the  memoryless  property  of  the  exponential  times  between 

course  changes .  The  results  of  such  a  simulation  are  pre- 

\  '  1 

sented  in  this  paper  in  the  form  of  improved  survivability 

for  the  evader. 

Five  course  change  rules  were  evaluate^  in  the  simula¬ 
tion.  The  probability  density  function  of  each  is  pre¬ 
sented  on  subsequent', pages  along  with  a  graph  of  the  re- 

1  \ 

i  suiting  attacker's  pk  versus  attacker  strength,  denoted  as 

S=A/ir(v»TL)  These  graphs  are  for  representative  values 

of  the  parameter  a.  Alpha,  a,  is  the  product  of  time-late 
\  and  the  exponential  parameter  X. 

•  -  —  - 

^Course  changes  were  assumed  to  be  independent  of  the 
underlying  Poisson  process . 
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The  simulation  established  a  new  lower  bound  on  the  pk 
the  evader  could  yield  to  the  attacker.  Figure  3  is  a  graph 
illustrating  the  improvement  of  the  new  bound  over  that 
found  by  Washburn.  This  new  bound  is  the  least  lower  bound 
of  all  the  versus  S  curves  for  the  four  values  of  a  and 
all  rules  simulated.  Examples  of  such  curves  are  shown  in 
Figs.  5-9.  The  old  bound  in  Fig.  3  is  the  lower  bound  of 
the  p^  curves  simulated  using  the  uniform  rule  exclusively. 

The  minimal  pk  as  a  function  of  ot ,  over  the  five  rules , 
is  shown  in  Fig.  4  for  three  specific  attacker  strength 
levels.  The  curves  in  Fig.  4  confirm  the  statement  in  Ref. 

2  that  the  evader  should  turn  most  often,  that  is  a  should 
be  highest,  against  the  weakest  opponent. 

Comparisons  of  the  p^  curves  of  the  different  rules 
have  shown  that  the  uniform  rule  is  not  optimal  for  all 
attacker  strengths  and  a  values.  For  example,  the  reverse 
course  rule  (see  Fig.  7)  was  shown  to  be  better  than  the 
uniform  rule  for  a*»l.  The  difference  between  these  two 
rules  is  illustrated  in  Fig.  10.  Against  a  weak  opponent 
(strength  <.6)  the  evader  would  do  better  to  use  the  re¬ 
verse  course  rule  rather  than  the  uniform  rule  for  o*l. 

This  improved  survivability  occurs  because  the  reverse 
course  rule  created  a  "flatter"  evader  position  density. 

None  of  'the  five  course  rules  evaluated  produced  a  con¬ 
sistently  smaller  pk  for  all  values  of  attacker  strengths 
and  a.  Hcwever,  the  left-right  rule  was  dominated  by  the 
other  four  in  all  cases .  Each  of  the  five  course  rules  was 


simulated  at  six  different  levels  of  a  by  varying  ET.  The 
graphs  of  pk  versus  A/ir(v*TL)2  for  each  simulation  are  pre¬ 
sented  in  the  Computer  Output  section  along  with  tabulated 
p^'s  for  specific  attacker  strength  levels  as  a  function 
of  the  course  rule  and  the  a  used. 


New  Bound  on  p^  That  the  Evader  Can  Yield 
as  a  Function  of  Attacker  Strength 

FIGURE  3 
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Density  Function  of  Modified  Left-Right  Rule 

FIGURE  6A 


Pk 


FIGURE  6B 
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\^x 

tt  2lt 

Density  Function  of  Reverse  Course  Rule 


FIGURE  7A 


Pk 


FIGURE  7B 
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Probability  {Course  Change  **  v/2]  -0.5 
Probability  [Course  Change  *  3ir/2]  -  0.5 

Density  Function  of  Left-Right  Rule 

FIGURE  8A 
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FIGURE  9B 
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a*»1.0  A/tt(vTL)2 

A  Comparison  of  the  Tradeoffs  Between  the  Uniform 
and  Reverse  Course  Rules  As  a  Function  of  Attacker  Strength 

FIGURE  10 
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IV.  SIMULATION  MODEL 


Let  E's  position  at  time  t  be  at  X=0,  Y*0  of  a  cartesian 
coordinate  system.  Then  rotate  the  axis  until  the  positive 
abscissa  is  aligned  on  and  pointing  in  E's  last  direction 
of  travel.  Then  regardless  of  the  course  rule  used  by  E 
there  is  at  least  exp (-a)  probability  that  he  will  be  at 
X«*v •  TL ,  Y«0  at  time  t+TL.  This  is  the  probability  that  E 
does  not  change  course  during  TL.  So  the  attacker's  prob¬ 
ability  of  killing  E  by  centering  a  small  portion  of  his 
lethal  area  at  (v*TL,0)  is  exp  (-a).  Then  add  to  exp (-a) 
the  integral  of  E's  density  function  over  that  part  of  the 
uncertainty  circle  where  the  rest  of  the  lethal  area  is 
targeted  to  determine  the  total  p^.  The  total  amount  of 
lethal  area  allocated  in  this  manner  is  the  numerator,  A, 
of  the  attacker's  strength  function  S. 

A  computer  program  to  simulate  the  above  computation 
was  written  in  FORTRAN  and  run  cn  the  IBM  360.  The  program 
consisted  of  two  phases,  first  the  play  of  a  strategy  and 
secondly  the  scoring  of  that  play. 

A.  SIMULATION  OF  A  STRATEGY 

To  create  a  single  play  the  evader's  track  was  simulated 
from  an  initial  position  at  time  t  to  the  resultant  position 
at  t+TL,  The  track  was  the  result  of  a  specific  maneuvering 
strategy  being  simulated.  A  strategy  was  made  up  of  two 
decision  rules.  The  first  rule  determined  the  times  be¬ 
tween  course  changes  which  were  exponential  random  variables 
with  mean  ET=1/X.  This  rule  was  common  throughout  all  of 
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the  simulations,  although  the  parameter  ET  was  a  variable. 
The  second  rule  determined  the  magnitude  of  subsequent 
course  changes.  The  only  other  kinematic  restriction  was 
a  constant  speed,  v,  for  the  evader. 

To  start  a  play  datum  was  initialized  by  setting  t=0r 
X(0)=0,  Y(0)=0  and  an  initial  course,  c^ ,  was  selected  from 
the  uniform  distribution.  Then  a  sequence  of  exponential 
times  [t^]^^  N(TL)^  were  generated  and  a  sequence  of  course 
changes,  n(TL)'  were  generated  using  the  course 

change  rule.  The  evader's  position  at  t=TL  was 


X(TL) 


N-l 

Z  V't^  (Cos(ci)-Cos(cN))  +  v*TL*Cos(cN) 


1.1 


Y  (TL) 


N-l 

E  v  t^  (Sin(ci)-Sin(cN) ) 


+  v*TL*Sin(cN) 


1.2 


B.  SCORING  THE  PLAY 

Once  E's  position  at  t=TL  was  determined  that  observa¬ 
tion  was  scored.  To  score  a  play  the  coordinates  of  E's 
position,  X (TL) ,  Y (TL) ,  had  to  be  transformed.  The  purpose 
of  the  transformation  was  to  make  the  observed  position  in¬ 
dependent  of  the  particular  initial  cource  c^.  The  trans¬ 
formation  was  a  rotation  of  the  coordinate  axis  about  the 
datum  so  the  positive  abscissa  would  be  aligned  in  the 
direction  cf  the  initial  course.  The  transformed  position 
was: 


6The  number  of  elements  of  this  sequence,  N (TL) ,  is  a 
Poisson  random  variable. 


22 


2.1 


X'  (TL)  =  X(TL)Cos{Cl}  +  Y(TL)Sin(c  ) 

Y'  (TL)  =  Y(TL)Cos(c1)  -  X(TL)Sin(c  )  2.2 

An  example  of  this  transformation  is  shown  in  Fig.  11. 

If  E  had  not  made  a  course  change  during  TL  then  his 
transformed  position  would  have  been 


X'  (TL)  *  vTL  3,1 

Y*  (TL)  =  0  3.2 

A  grid  system  of  square  cells7  was  placed  over  the 
playing  area  of  Fig.  11B  and  a  determination  was  made  as  to 


ORIGIN;  Datum 

B;  X (TL) ,  Y (TL) 

C^;  Initial  Course 

SAMPLE  PLAY  PRIOR 
TO  TRANSFORMATION 


ORIGIN:  Datum 

B:  X' (TL) ,  Y' (TL) 

A:  (v-TL,0) 

SAMPLE  PLAY  AFTER 
THE  TRANSFORMATION 


FIGURE  11A  FIGURE  11B 

which  cell  X* (TL) ,  Y 1 (TL)  was  in.  Each  cell  of  the  grid 
had  an  associated  value  which  would  represent  the  number  of 
times  a  play  resulted  in  an  observation  in  that  cell.  When 


7 

The  cumulative  cell  area  was  normalized  by  the  factor 
TT  (  V*  TL)  c  . 
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I 


the  appropriate  cell  wag  determined,  for  the  play  being 

V  I 

scored,  that  value  was  incremented  by  one. 

A  simulation  run  was  composed  of  16810  plays  and  scoring 
iterations  for  a  strategy  utilizing  a  specific  course  change 
rule.  The  input  variables  for  a  run,  besides  the  course 

| 

rule  wore  v,  TL  and  ET.  At  the  end  of  a  simulation  run  the 
grid  system  was  a  two  ^imensional  histogram  of  E's  position. 
The  frequencies  in  the  histogram  were  then  ordered,  accumu¬ 
lated  and  normalized  to  achieve  cumulative  cell  probabilities, 
The  ordering  corresponded  to  the  conservative  assumption  that 

P  could  divide  hie  lethal  area  and  target  only  those  cells 

\ 

with  the  higher  probabilities.  These  probabilities  were 


then  plotted  agAinst 
8 


the  cumlative  cell  area  they  repre¬ 


sented.  This  graph  was  labeled  p^  versus  attacker  strength^ 

and  was  the  primary  output  of  the  program. 

\  ^ 

Thirty  simulation  runs  were  made  to  investigate  five 
'different  course  change  rules  and  six  different  values  of 
ET.  For  all  luns  the  following  constant  values  were  main¬ 
tained,  v=5  knots  and  TL=2  hours.  The  graphs  from  those 

\ 

runs  Are  included  in  the  Computer  Output  section.  Each 
course  change  rule  was  simulated  many  different  times,  the 
only  difference  between  runs  being  the  pseudo  random  num¬ 
ber  generator  seeds.  This  was  done  to  check  for  the  vari¬ 
ability  of  the  pfc  graph  for  that  rule.  In  all  such  runs 


The  amount  of  area  in  a  cell  was  denoted  as  cell  size  and 
equal  to  '  (2»TL/41) 2.  TL  and  v  were  held  constant  there¬ 
fore  cell  si'-se  was  always  approximately  1/4  nm2. 
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the  resultant  graphs  were  so  similar  that  any  difference 

l 

was  not  distinguishable.  For  simulation  run  required  two 
minutes  and  forty- five  {seconds  using  125  K  on  the  IBM  360. 

| 

\ 

1 

I 

t 


V.  CONCLUSION 


It  was  determined  that  the  uniform  course  change  rule 
is  not  optimal  within  the  Poisson  class  of  strategies.  Com¬ 
parison  of  various  rules,  see  Fig.  10,  showed  that  under 
certain  conditions  a  single  rule  such  as  the  reverse  course 
change  rule  is  better.  Also  the  new  lower  bound  on  p^, 
taken  over  all  rules  and  values  of  a  simulated,  is  an  im¬ 
provement  on  the  lower  bound  achieved  from  only  the  uniform 
rule. 

Certainly  not  all  the  possible  course  change  rules  were 
simulated.  The  rules  evaluated,  however,  were  representa¬ 
tive  of  the  broad  class  of  possible  rules.  When  compared 
to  the  uniform  rule,  all  but  one  cf  uie  other  rules  showed 
that  E  could  improve  his  situation  if  he  knew  P's  strength 
by  selecting  the  better  rule  for  that  encounter.  The  one 
rule  that  was  consistently  dominated  was  the  left-right 
rule. 

It  remains  unknown  whether  or  not  the  optimal  strategy 
is  a  Poisson  strategy. 
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APPENDIX  A 


FLOW  DIAGRAM 


l  nil  . . . .  ■ 

- ipseudo-random  number  seed,  an  input  variable 


- i evader  speed  in  knots 

- itime-late  in  hours 


- {mean” time  between  course  changes,  an  input 

{variable 

- {playing  area  and  histogram;  each  cell  is 

ja  square  with  dimension  (2*V*TL)/41 


- i INDEX  is  a  book-keeping  record  of  PP  cell 


|f easibility 


PP(i,j)  an  exterior  cell  \ 

_ Sg _ 

PP(i,j)  an  interior  cell  | 


No 


O 
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Generate  a  course 
change  c  by  the 
rule  being  sim- 
lated 


30 


33 


Rule 

1.0 

COMPUTER  OUTPUT 

TL/ET 

2.0 

3.0 

4.0 

Uniform 

<729 

.580 

.501 

.486 

Modified 

Left-Right 

.737 

.580 

.501 

.486 

Reverse 

Course 

.690 

.572 

.548 

.564 

Left-Right 

.917 

.768 

.650 

.588 

Truncated 

Uniform 

.697 

.556 

.509 

.517 

column 

minimum 

.690 

.'556 

.501 

.486 

for  attacker  strength  equal  0.2 
Table  I 
TL/ET 


Rule 

1.0 

2.0 

3.0 

4.0 

Uniform 

.854 

.784 

.752 

.760 

Modified 

Left-Rlaht 

.878 

.  7  84 

.752 

.768 

Reverse 

Course 

.838 

.791 

.799 

.831 

Left-Right 

.964 

.909 

.870 

.854 

Truncated 

Uniform 

*.846 

.768 

.768 

.815 

column 

.838 

.768 

.752 

.760 

minimum 

Pk  for  attacker  strength  equal  0.4 
Table  II 
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Course  Change  Rule:  Uniform 


0.2  0.4  0.6  0.8  1.0 


A/xCv-TL)2 

Pk  versus  Attacker  Strength 

Course  Change  Rule:  Uniform  El  *>  0.66 


0.2  0.4  0.6  0.8  1.0 

A/TV(v.TL)2 

Pk  versus  Attacker  Strength 

Course  Change  Rule:  Uniform  ET  *  0,5 


42 


0.2 


0.4 


0.6 


0.8 


1.0 


Pk  versus  Attacker  Strength 
Course  Change  Rule:  Modified  Left-Right 


A/*  (v.TL)2 

ET  -  2.0 


43 


0.2 


0.4 


0.6 


0.8 


1.0 


A/lt(v-TL)2 

Pk  versus  Attacker  Strength 

Course  Change  Rule:  Modified  Left-Right  ET  *■  0.66 


versus  Attacker  Strength 
Course  Change  Rule:  Modified  Left-Right 


0.2  0.4  0.6  0.8  1.0 

A/H  (v.TL)2 

versus  Attacker  Strength 

Course  Change  Rules  Reverse  Course  ET  ■  1.33 


i 


0.2  0.4  0.6  0.8  1.0 

A/TC  (v*TL)~ 

Pk  versus  Attacker  Strength 

Course  Change  Rule:  Reverse  Course  ET  »  1.0 
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Pk  versus  Attacker  Strength 

\ 

Course1 Change  Rule:  Reverse  Course  ET  ■  0.8 


\ 


\ 

I 

i 


52 


A/lf  (vTi.) 


\ 


versus  Attacker  Strength 

i 

Course  Change  Rule:  Reverse  Course  \  ET  «  '0.66 

\  : 


0.2 


0.4 


0.6 


0.8 


1.0 


A/n  (v*TL)  2 

versus  Attacker  Strength 

Course  Change  Rule:  Left-Right  ET  «  2.0 


0.2 


0.4 


0.6 


0.8 


1.0 


A/H(v.TL)2 


Pit  versus  Attacker  Strength 

Course  Change  Rule:  Left-Right  ET  ■  1.33 
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0.2 


0.4 


0.6 


0.8 


1. 


A/  7t  (v*  TL)  2 


versus  Attacker  Strength 
Course  Change  Rule:  Left-Right 


ET  -1.0 
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0.2 


0.4 


0.6 


0.8 


versus  Attacker  Strength 
Course  Change  Rules  Left-Right 


A/H  (v.TL)2 


ET  -  0.8 
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versus  Attacker  Strength 
Course  Change  Rule:  Left-Right 


0.2 


0.4 


0.6 


0.8 


A/n  (vTL)2 


versus  Attacker  Strength 

Course  Change  Rules  Truncated  Uniform  F.T  «  2.0 
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0.2 


0.4 


0.6 


0.8 


1.0 


A/H  (vTL)2 


versus  Attacker  Strength 

Course  Change  Rule:  Truncated  Uniform  ET  ■  1.33 


62 


A/It  (v.TL)2 

Pk  versus  Attacker  Strength 

Course  Change  Rule:  Truncated  Uniform  ET  »  0.8 


64 


A/tt  (vTL) 


versus  Attacker  Strength 
Course  Change  Rule:  Truncated  Uniform 


ET  «  0.66 
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Course  Change  Rule:  Truncated  Uniform  j  ET  «  0.5 

\ 

I 
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uouu  ouuoouuouDyuuuuuowu  uuoouu  uo  uu  oooc;;yy 


COMPUTER  PROGRAM 


\ 


A  SITUATION  OF  THE  CONTINUOUS  EVASIVE 
GAME  FOR  THE  POISSON  CLASS  OF  STRATEGIES 


INTEGER 

INTEGER 

INTEGER 

INTEGER 

INTEGER 

INTEGER 


*  2  PP.\lNDEX 

*  2  .CUMSUM 

*  2  NCQUNT 
*2  XP.QQ 

*  4  C.RITE3, 

*  4  TOTCEL 


\ 


REAL  *  4  IN 
REAL  LABEL/'  '/ 

REAL*a  ITITLE{12)/* CLOTHIER* • * 

X  *,' STRATEGY' ,  •  TRNCRSS 

X2*'  4.0'/ 

BY  MAKING  THE  HISTOGRAM  INTEGER  STAR  TWO 
OF  STORAGE  REQUIRED  FOL  THE  PROGRAM  WILL 


THOMAS'.'  J, 

'  (LAMDAS  •  * 


'  ,3*» 

TL  )-  % 


THE  AMOUNT 
BE  REDUCED. 


DIMENSION  PP(41,41  I  .INDEX ( il ,41 ) 

DIMENSION  XP(l68ii ,00(41,41) 

EQUIVALENCE  < PP ( 1. 1 ) , XP( I ) ) 

0|  MANSION  EEX(2) , 5 =Y ( 2 ) , EX < 692 ) , EY ( 692 ) , 6XX ( 692 ) 


;NSI  ON  E  YY(  692) 


\ 


PP  IS  THE  PLAYING  AREA  AND  HISTOGRAM.  EACH  CELL  IN 
PP  REPRESENTS  AM  AREA  OF  SIZE  C5  LSI  Z.  INDEX  !S  A 
BOOK-KEEPING  ARRAY  USED  T0\ DETERMINE  WHICH  CELLS  IN 
PP  ARE  FEASIBLE,  E.G.  IF  INDEX(I,.J)=0  THEN  PP  <  I ,  J  )  IS 
AN  INFEASIBLE  CELL  •  SIMILARLY  IF  INDEX(  I  ,  J)=l  OR  2 
THEN  PP(1,J)  .IS  A  FEASIBLE  CELL. 


PP  IS  A  SQUARE  WITH  THE  TIME-LATE  CIRCLE  INSCRIBED  IN 
!T,,  THEREFORE  pP  HAS  SOME  CELLS  THAT  ARE  NOT  INSIDE 
THE  TI MS-LATE  CIRCLE  AND  FOR  THOSE  CELLS  THEIR  INDEX 
VALUE  IS  ZERO.  FOR  THOSE  CELLS  IN  PP  THAT  ARE  CON¬ 
TAINED  ENTIRELY  IN  THE  CIRCLE  THE  ASSOCIATED  INDEX 

VALUE  IS  ONE.  * -  -  ~~  “ 

THAT  ARE  ON  THE 
EDGE  CELLS  AMD 


i  IN  THE  CIRCLE  THE  ASSOCIATED  INDEX  1 
AND  LASTLY  THERE  ARE  THE  C=LLS  IN  P° 

E  EDGE  Of-  THE  CIRCLE,  THESE  ARE  CALLED 
HAVE  THE  INDEX  VALUE  OF  TWO. 


NCCUNT=0 
PI=3i 1415926 


I 


DATA  iXjVjTLtET/STDlStOB.OfOa.OOtOO.B/ 


i  A  ,  V  J  ILlI 

V  IS  PARTI  CLG  SPEED 
TL  IS  THE  TIME-LATE 

ST  IS  THE  MEAN  TIME  BETWEEN  COURSE  CHANGES 

IX  IS  THE  SEED  FOR  THE  PSEUDO-RANDOM  NUMBER  GENERATORS 


\  DSLTAX=(  2.  0*V*TL  1/41.0 
1  DELTAX  IS  THE  HORIZONTAL  DIMENSION  OF 


i 


A  CELL. 


CELSI Z-OELTAX**  2 

CELS  IZ  IS  THE  AMOUNJ  OF  REAL  AREA  CONTAINED  IN  A  CELL. 


K  =  ?0 

THE  NUMBER  OF  CELLS  IN  EITHER  DIMENSION  IS  ALWAYS  A 
CONSTANT  41.  THIS  CONSTANT  IS  EXPRESSED  AS  AN  ODD 
INTEGER  AND'i  WRITTEN  IN  THE  FOLLOWING  MANNER, 

41=  2  *  <K)  +  1.  THEREFORE  K  =  20. 


ZERO  OUT  THE 
00  5  1=1,41 

DC  4  J  =  1 *  41 
INDEX! I, J )=0 


HI STOGRAM 
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OUUOUUU  UUOU  UJO  OUO  UOUOUUOOUOU 


PP(I,J)=0 

QQ(I,J)*0 

4  CONTINUE 

5  CONTINU* 


NOW  TO  CHECK  EACH  OP  THE  CELLS  IN  THE 
DETERMINE  ITS  ASSOCIATED  INDEX  VALUE. 


HISTOGRAM  AND 


41 

42 


QUADRANT  I 
DO  42  J=21 *41 
DO  41  1=1,21 

0UT=(V*TL-J*DELTAX)**2+(  (I  -1 1  * DELTAX-V4T U **2 
IN  =(V*TL-<  J-l  >*D5LTAX)**2+(I*DELTAX-V*TL) **2 
RADSQ= ( V*TL )**2 

I F I C  OUT  «LE  «  PAD  SOI.  AND.  (IN.  LE. RADSQI  JINDEXtI.J  1=1 
IF  <  (OUT.GT.PADSO).  AND  .  ( IN.  GT.R  ADSQ  ))INDEX(  I,  J)=0 
I  PC ( OUT. GT .RADSQi • AND. (IN. Lc.RADSO ) ) INDEX ( I, J)  -2 
CONTINUE 
CONTINUE 


43 

44 


QUADRANT  II  CHECK 
DO  44  J=] ,21 
DO  43  1=1.21 

OUT=(V*TL-<  1-1  )*DELTAX)**2+(V*TL-(  J-l >*DEL TAX  1**2 
IN=(  V*TL-I*D?-L TAXI **2+(V*7L-J* D5LTAX  1**2 
RADSO=(V*TL 1**2 

I F  C ( OUT. LE. RADSQI. AND. ( I N. LE .RADSO I )INDEX( I,J)  =  1 
I F( ( OUT .GT.R AD  SO  I  •  AND. ( I N. GT. RAD  SO ) ) I ND2  X ( I » J ) =0 
IF ( (OUT . GT .RADSO I . AND. CIN.LE.RADSQ  I  )INDEX( I, J 1  =  2 
CONTINUE 
CONTINUE 


45 

46 


QUADRANT  III 
DO  46  J=l. 21 
DO  45  1=21,41 

OUT=(V*TL-(  J-l)*DELTAX)**2  +  <I*DEl.TAX-V*TL)**2 
IN=(V*TL-JTDELTAX)**2+( (I-1)*DELTAX-V*TL )+*2 
RADSQ=( V*TL)**2 

IfTcoUT.LE.RADSQ  I.  AND.(  IN.LE.RADSOI  )IN0£X(  I .  J)=l 
I F( (OUT. GT.R ADSQ I. AND. ( I N.GT. RADSO)  nND6X(  If  J  )  =  0 
IF(  ( OUT. GT.  RADSQI.  AND.  (I  N.  LE.  RADSO  I  )  I  NDSX  (  I.  j)=2 
CONTINUE 
CONTINUE 


47 

48 


QUADRANT  IV  CHECK 
DO  48  J=21 *41 
DO  47  1=21.41 

OUT*(  VtfTL-l'^DE LT AX  )**2+ ( V*TL-J*DELTAX 1**2 

IN  =(  V*TL-(  l-l)*DELTAXH*2+{  V*TL-(  J-l  I* DEL TAX  1**2 

RADSQ= ( V*TL  >*T? 

IP( (OUT. LS. RADSO I. AND. (IN. LE. RADSQI ) INDEX { I, J 1=1 
IF((  OUT,  GT.R  ADSQ).  AND.  (  IN.  GT.RADSQ)  )IND£X(  I  .J)=0 
IF ( (OUT. GT. RADSQI. AND. (IN. LE. RADSO I) INDEXt I, J 1=2 
CONTINUE 
CONTINUE 


NOW  TO  SIMULATE  THE  PARTICLE  MOTION  UNDER  THE 
SPECIFIED  RULES.  THIS  IS  THE  MAIN  DO  LOOP  IN  THE 
PROGRAM.  FIRST  A  TIME  UNTIL  THE  NEXT  COURSE  CHANGE 
IS  GENERATED.  THEN  A  COURSE  TO  BE  STEERRED  IS  GEN¬ 
ERATED.  XX  AND  YY  ARE  THE  TWO  COMPONENTS  OF  POSITION 
ADDED  DUE  TO  A  COURSE  AND  TIME  SEGMENT.  X  AND  Y  AR 8 
THE  UPDATED  POSITION  cROM  DATUM.  AT  TIME-LATE  THE 
UPDATED  POSITION  IS  CONVERTED  TO  A  CELL  POSITION  AND 
THAT  CELL  VALUE  IS  INCREMENTED  BY  ONE. 
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c 

c 


DO  3  J»l. 16810  < 

X=0.0 

Y=0.0 

CLOCK*0.0 

C=0.0 

TR=TL-CLOCK 

CALL  RANDUUX,IY,R) 

IX*IY 

CALL  CORSET R »C ) 

THETA *C 
GO  TO  66 
1  TR=TL-CLOCK 

CALL  RANDU<!X,IY,R1 
IX*IY 

CALL  TRMCRS  (R  *C) 

...66 -CALI  RANOUU  X,  I  Y.,R  I 
IX  =  IY 

CALL  TRKTM<?T,R,T) 

TT=AMIM1 ( T  »TR) 

CALL  MOV«(V,TT,C,XX,YY) 

CLOCK=CLOCK+T 

X=X+XX 

Y=Y+YY 

IF(CLOCK.LT.TL)GO  TO  1 

XPRI ME=  X*CO  S ( THE  TA ) +Y*SIN< THETA  I 

YPRI ME=Y*COS( THETA ) -X^S IN (THETA ) 

X=XPRI ME 

Y=YPRIM= 

CALL  RANOIM  IX.  IY  »R  ) 

I  X=I  Y 

CALL_HisTO(K*  DELTAX  » R . CORK  *  X . i > 

CALL  RANPUUX.  IY.R) 

IX  =  IY 
Cl1RRs“  1  0 

CALL~HISTO(K, DELTAX. R,CORR,Y,I ) 

THE  PARTICULAR  CELL  I.J  IS  DETERMINE 
IV£  CALLS  OF  HISTO.  IT  IS  POSSIBLE 
COULD  BE  DETERMINED  TO  Be  IN  AN  INFS 
IF  THIS  HAPPENS  THEM  RESOLV  WILL  BE 
THAT  CELL  ASSIGNMENT  IN  A  PRESCRIBED 
FEASIBLE  CELL.  _  ^ 

I F ( INDEX (M.L ) .EQ.O) GO  TO  31 

30  PP(M,U  =  PP(M,L  )+l 
GO  TO  3 

31  CALL  RESOLV<L.M, I.N) 

NCOUNT  =NC0UNT+1 

L=I 

M=N 

GO  TO  30 
3  CONTINUE 

NUPCEL=0 
NUPEDC=0 
CRITS2=0.0 
DO  52  1=1 .41 
DO  51  J=l»41 
I F ( PP (  I , J) .EQ.OJGO 
NUPCEL=NUPCcL+i 
I F {  INDEX! I . J ) .NE • 2) GO  TO  51 
NUPEDC=NUPSDC+ 1 
CRITE2=CRITC2+PP(I »J)/1681>  .0 

51  CONTINUE 

52  CONTINUE 
CRI  TE3  =  NUP’fDC 
RNUPED=N UPSDC 
RNUPCE=NUPrFL 
CRITE1=RNUPED/RNUPCE 


TO  51 


D  BY  TWO  SUCCESS- 
A  POSITION  X i Y 
ASIBLE  CELL.  „ 
CALLED  TO  CHANGE 
MANNER  TO  A 
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fit 


204 

205 

206 


207 

206 

209 


C 

C 

C 

C 

C 


210 

211 


NUP1NC=NUPCEL-NUPEDC 
VIR1TF  <  6  *  203 1 

FORMAT ( ' 1 1 ,60X , *  PP  PAGE  1 • ) 

00  205  1=1.41 

HR!TE(6,204)  (  PP ( I  * J ) , J* l » 15 ) 

FORMAT  ( •  '  ,1517) 

CONTINUE  ^  4 

FORmIt( '{^60X,*PP  PAGE  2') 

DO  208  1*1.41 

WRITE (6, 207)  (PPM ,  J ) ,  J*15 ,29 ) 

FORMAT(*  *.1517) 

CONTINUE 
WRITF  ( 6,209) 

FORMAT ( • 1 • . 60X. 'PP  PAGE  3') 

DO  211  1=1,41 

WRITE(6»210)(PP(I,J)»J=29,41) 

FORMAT  (•  ',..131 7) 

CONTINUE 

NOW  THAT  THE  HISTOGRAM  IS  BUILT  ITS  ELEMENTS  ARE 
ORDERED  SO  THE  HIGHEST  VALUES  ARE  IN  THE  LOWEST  NUM¬ 
BERED  CELLS* 


XPd+l n  GO  TO  1010 


S 

c 

c 


1000  IC=0 

DO  1010  1=1.1680 
iFtxpm.GE. — 

MAX=XP ( 1+1) 

XP(I+1I=XP<I) 

XP(I)=MAX 

IC=10 

1010  CONTINUE 

I F ( I C. NE . 0 ) GO  TO  1000 
DO  2000  1=1. 4.1 
DO  2010  J  =  1 » 41 
QQ ( I ,J) =PP( J.I ) 

CONTINUE 
CONTINUE 
DO  1030  1=1,41 
DO  1020  J= 1 « 41 
PP(I. J)»00(I,J  ) 

CONTINUE 
CONTINUE 

AFTER  ORDERING  THE  CELL  VALUES  THEY  ARE  ACCUMLATGD  SO 
EACH  CELL  CONTAINS  THE  SUM  OF  ITS  OWN  VALUE  AND  ALL 
THE  VALUES  PRECEDING  IT. 


20 10 
2000 


1020 

1030 


c 

c 


10 

11 


32 

33 


CUMSUM=0 

DO  11  1=1,41 

DO  10  J-1,41 

CUMSUM=CUMSUM+PPn  ,J) 

PP(I , J)=CUMSUM 

CONTINUE 

CONTINUE 


NUMEXC=0 
NUMI NC  =0 
NUMSDC=0 
DO  33  1=1,41 
DO  32  J= 1 , 41 

IF  (INDSXi I, J1 *EQ.O)NUMgXC«NUMSXC+l 
IF  (INDcXM  ,.l)  .EO.l )  MJMI NC  =  NUM INC+1 
IF  (  INDEX(I.J».E0.2)NUMcDC=NUMEDC+l 
CONTINUE 
CONTINUE 

NUMCEL=NUMINC+NUMEDC 
TOTCEL  =  NUMEXC+NUM£  DC+NUMINC 
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oo 


NOW  TO  LOAD  THE  FIRST  900  CELL  VALUES  INTO  THE  PLOT¬ 
TING  VECTOR  EY.  THE  CUMULATIVE  NUMBER  OF  LOOKS  OR 
KILL  POWER  WILL  BE  PLOTTED  ON  THE  X  AXIS. 


■X(1)=0.0 
=X  C  2 1*1-0 
|Y( 1) =0. 0 
:Y( 21=  1.0 


DO  15  1=1,16 
DO  14  J=l, 41 
N={f-l)*4i+J 

fY(N)=PP(I*J)/1681G.O 
X1N)»N/ 1-3  85*0 
CONTINUE 
CONTINUE 
DO  16  J=1 *36 
1*17 

N=(I-1)*41+J 
EY(NI=PP(I ,J)/16810.0 
EXIN )=N/13  85.0 
CONTINUE 


DO  17  J*37.41 
1=17 

N=(I-1)*41+J 

L=N-692 

EYY(L)=PP( I, JJ/16810.0 
E  XX( L ) =N /1385. 0 
'  CONTINUE 
DO  19  1=18,33 
DO  18  J=1.41 
N=(T-1)*41+J 
L=N-692 

EYY(U  =  PP(  I.JJ/16910.0 
EXX(L)=N/1385.0 
}  CONTINUE 
)  CONTINUE 
DO  20  J=1 ,31 
1=34 

N=( I-l)*41+J 
L=N-692 

EYY(L)=PP( I • J ) /16810.0 
EXX ( L ) =N/1385.0 

>  CONTINUE 

WRITE(6,606)  NCOUNT 

>  FORMAT!//, 25X, 'NCOUNT  EQUALS  ',116) 

WRITE  (6,100)  NUME XC , NUMEDC , NUMI NC ,TOTCEL 

)  FORMAT  (/.10X,  'NUM£XC=«,  116,  5X,  'NUMEDC*'  ,  II  6, 5X, '  NUMIN 
XC=' ,11 6,5X,'T0TC=L=' ,116) 

WRITE! 6, 101 )  CRITS3,CRITE1,CRITE2 
L  FORMAT  ( //, 10X , ' NUMBER  OF  POSITIVE  EDGE  CELLS  IS  ',11 
X6,//,10X,'THE  PERCENTAGE  OF  POSITIVE  CELLS  WHICH  AR*  E 
XDGF  CELLS  IS  • , IF  4. 3, // , 10X , • THE  AMOUNT  OF  PRQQABI LIT 
XY  IN  THE  EDGE  CELLS  IS  *,1F4.3) 

CALL  DRAW  (  2  ,EE  X ,EFY  ,1  ,0  , LABEL ,  ITITLE ,0  ,0  ,0 ,0  ,0 , 0, 9,  9, 
XO.LAST ) 


CALL  DRAW 
XO.LAST) 
CALL  DRAW 
X9|0jLAST) 


;uvjc  io  !  irtij) 

(  2,EEX,EFY,1  ,0,  LABEL,  I  TITLE, 0,0, 0,0  ,0,0, 9,  9, 
(692  ,  EX,  £Y,  2,0,  LABEL,  ITITLE,  0,0,  0,0,  0,0,  9, 9, 
(692, EXX, EVY, 3,0, LABEL, I  TITLE ,0,0, 0,0, 0,0, 9, 
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SUBROUTINE  RANDUU  X,IY.YFL) 

1Y=IX*65539 

IF(IY) 5,6.6 

5  I  Y  =  I Y+2 1.47483647+1 

6  YFL**IY 

YFL“YFL*.4656613E-9 

RETURN 

END 


SUBROUTINE  TRKTM  (E  T  .R .  T) 

RA=1.0-R 

T1=-ET*AL0G(R) 

T2=-ET*AL0G(RA) 

T=T1 

RE-TURN 

END 

SUBROUTI N«  MOVE ( VS L ♦ TI ME , ANG , X  , Y ) 
X=VEL*TIMg*COS(ANG) 

Y=VEL*TI ME*SIN (ANG ) 

RETURN 

END 

SUBROUTINE  HI STO (K , OELTAX, R,COP.R tX , I ) 

UPBND= ( ( 2*K+1 )*OEL TAX ) /2. 5 

L=2*K+1 

XX= ( X+ ( CORR+UPBNOJ ) /OELTAX 

XX=CORRfXX 

XI  =AMOD  (  XX  ,1  .0  l 

M=INT( XX ) 

IF ( M . EQ.O I  GO  TO  1 
I F (  M.EO.UGO  TO  2 
I F (  X 1  .NE  .O.OGO  TO  3 
IF ( R. GT  .0  •  5  I  GO  TO  3 
I  -M 

GO  TO  4 
1  1  =  1 
GO  TO  4 

2  I=L 

GO  TO  4 

3  1  =M+1 

4  RETURN 
END 


SUBROUTINE  RE SOIVI I , J . L . MJ 

IFII.GT.2UG0  TO  2 

IFtJ.GT.2UG0  TO  1 

L=  I  +  1 

P=J  +  l 

GO  TO  4 

L=I  +  1 

M=J-l 

GO  TO  4 

IF( J ,GT .21 )G0— TO  3 
1  =  1-1 
P=J+1 
GO  TO  4 
1  =  1-1 
M  =  J- 1 
RETURN 
END 
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SUBROUTINE  CORSE(R,CJ 

PI»3. 1415926 

A*0.0 

B=2.0*P! 

F«2«0*PI 

DELTAC«R*(B-A)+A 

C*£+DcLTAC 

C*AMOD<C»F) 


c 

c 

c 

c 

c 

c 


c 

c 


c 

c 


c 

c 


ruis°theNuniform “variate  used  to  generate  the  next 
COURSE  CHANGE. 

IS  THE  PRESENT  COURSE  AND  WILL  BE  RETURNED  AFTER 
DOING  £  NEW  HEADING  CHANGE ■ 

DELTAC  IS  THE  AMOUNT  OF  THE  NEW  HEADING  CHANGE 


X 


S=0»  20 


S  IS  THE  TEEPEE  SLOPE  AND  IS 
2/<PI**2)  =  0.202.  NOTE  ,  S 


BETWEEN  0  AND 
CANNOT  EQUAL  0. 


PI =3. 141 5927 

IF( 0.0. LE.R .AND. R.LE.0.25)G0  TO  1 
IFI0.25.LT. R. AND. R. LS.0.5)GO  TO  2 
IF(0.5.lT.R.ANO.R.LE.O.75)G0  TO  3 
IF(0.75. LT.R.ANO.R. LE.l.OJGO  TO  4 


1  A=S/2.0 

B=<1.0/{2.0*PI ) )-<  S*P I /4.0) 

D=-1.0*'-R 

E  =  B**2-4.0*A*D 

PELT  AC= ( +SQRT (Etl/I 2.0*A) 
GO  TO  5 


RPRI ME=R— 0 .5 
A=  S/2. 0 

8=‘1.0/( 2.0  +  PI  )  )-<  S*PI /4.0) 

D=-l .0*RPRIM£ 

E=B**2-4.0*A*0 

DELT  AC= ( -1 .0*8  +SQRT { E  ) )/( 2.0*A) 

DELTAC=DELTAC+PI 

GO  TO  5 


2  A=(-1.0*S)/2.0 

B=(1.0/(  2.0*PI  ) )  +(  3,0*S*PI/(  4.0)  ) 
D=“l  .0*(  ( S*P  I*  2  1/4,0  )  -R 
E=B**2-4.0*A*D 

DELTAC=<-1.0>*B+SQRT<S  ))/(  2.0*A) 

GO  TO  5 


4 


RPRI  MS=R“‘0 .5 
A“(-l.  0*  S)  /2.0 

B~  ( 1 .0/ (  2.D*P  I  ))+(  3.0*S*PI/<  4. OH 
0*-l  ,0*<  (S*PI**2  )/4.0  ) -RPRI ME 
E=B*72-4.0*A*D 
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DELTAC»(-1*0*B+SQRT <E) )/(2.0*A) 
OgLTAC*OELTAC+PI 

5  C-C+DEITAC 

F=2.0*PI 

C*AMOD(C » F) 

RETURN 

END 


SUBROUTINE  PC0RS5(R,C> 

PI  =3. 1415926 
F=2.0*PI 
,  S«1 *0/ ( P!**  > 

IFIR.GT.O.  B)GO  TO  1 
A=2.0*R/S 
OELTaC=SQRT ( A) 

GO  TO  3 
A=-1.0*S/2.0 

8=(  1.0/  (  2»0**PI  .  +(3.0*PI*S/2.0) 

D»-1.0^S*S(P!*«  M-R 

E»(8**2I-4#0*A*D 

0ELTC1  =  (-B+SCjRT(  =)  )/<2.0*At 

DELTC2=(-B-S0P.T(  =  n  /<2.0*A) 

IP(PI .LT.DELTC1.AM0.DEITCI.LE.F)G0  TO 

06LTAC=DELTC2 

GO  TO  3 

DELTAODSLTC1 

C=C+D2LTAC 

C=AM0D(C,F) 

RETURN 

END 


+<3.0*PI*S/2.0) 
)  -R 


SUBROUTINE  LRCORSIR.CJ 
PI=3. 1415926 
F-2»0*PI 

1FIR.GT.0.5IG0  TO  l 
C=C+ { PI /2 .0  ) 

GO  TO  2 

C=C+(  <3.0^1  1/2,0 » 
C=ANC1D{C»F) 

RETURN 

END 


SUBROUTINE  TRNCRS(R,C) 

PI= 3. 141 5926 

F=2.0*PI 

A=Pl/6.0 

B= ( 11 *0  *PI J/fe.O 

0ELTAC=R*(8-A»  *A 

C=C+DELTAC 

C=  AMOD(  C »  F  ) 

RETURN 

END 
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